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1974 Annual Report on the Investigation 
of Critical Burning of Fuel Droplets 


Summary 


This report discusses activities under NASA Contract NGR 39-009-077 for 
the period January 1, 1974 to December 31, 1974# During this period the work 
was divided into two phases: 1. High pressure droplet burning in flowing 

combustion gases and 2. Combustion response of monopropellant droplets. The 
results for each phase of the study for the present report period may be 
summarized as follows: 

1. High Pressure Droplet Burning in Flowing Combustion Gases . Theory and 
experiment were compared over the currently available range of test data 
for droplet combustion in combustion gas environments. Predictions were 
within ±30% of measurements for a wide range of conditions when the model 
employed a multiplicative correction for convection based on average 
properties between the liquid surface and the flame. The range of 
comparison for these results is summarized in Table 1. With regard to 
liquid surface conditions, the low pressure model is adequate at pressures 
less than one half the pressure required for critical combustion, while 
the high pressure version provides a superior prediction at higher pressures. 
The two models do not overlap, however, within their range of validity. 
Therefore an intermediate pressure model f which is being developed to provide 
greater accuracy in this region, is presented in this report. 


2. Combustion response of monopropellant droplets . Earlier response results 
for burning liquid monopropellant strands were applied directly to the 
discussion of the response of large droplets. The limits for the ’’large 
droplet” conditions are discussed. The results show sufficient response 
to provide a viable mechanism for combustion instability; in the 
frequency range associated with combustion instability and for droplet 
sizes present in sprays. The results also suggest that gas phase transient 
phenomena might contribute sufficient gain for combustion instability at 
low pressures. Thus although the gas phase might be considered quasisteady 
for lifetime calculations, important contributions to the response could be 
overlooked by this approach. The strand analysis has been converted to 
consider the combustion response of droplets under the same property 
assumptions. This model allows evaluation of the limits of the ’’large 
droplet” model to be determined as well as combustion response down to the 
limit of low reactivity (evaporative limit). Zero-order closed form 
solutions have been obtained for the evaporative limit. Numerical calcu- 
lations are being conducted for the response and zero-order results 
for finite reactivity. A quasisteady gas phase model has also been 
completed. 
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1974 Annual Report on the Investigation 
of Critical Burning of Fuel Droplets 


1. Introduction 

The objective of this study was to continue earlier work on the steady and 
unsteady combustion of liquid fuel droplets under rocket engine conditions. 

Emphasis has been placed on consideration of combustion at elevated pressures 
and temperatures representative of typical combustion chamber environments. 

In addition, the problem of the open- loop response of burning fuel droplets to 
imposed pressure oscillations is being investigated. The results of this study 
has technical application to the design of liquid fueled rocket engines and 
the determination of the combustion instability characteristics of these 
engines. 

This report gives a summary of progress on the investigation for the period 
January 1, 1974 to December 31, 1974. During this period, the work was divided 
into two major phases, as follows: 

1. High pressure droplet burning in flowing combustion gases. Earlier work 
provided data on droplet gasification rates in flowing combustion gases 
over a wide range of temperatures, pressures, ambient oxygen concentra- 
tions, fuels, and Reynolds numbers (1,2). This data was analyzed in 

an attempt to improve the prediction of gasification rates from the 
results of Ref. 1. Improvements have been made by modifying the 
correlation for convection effects. A second phase of this portion of 
the study involved developing an intermediate pressure, phase equilibrium 
model in order to obtain a continuous transition between the present 
low and high pressure versions. 

2. Combustion response of monopropellant droplets. This portion of the 
study is concerned with extending the analysis of the "open loop response 
of burning monopropellant strands (3), to the case of droplets. The 
basic method employs a perturbation analysis allowing for both liquid and 
gas phase transient effects. In addition, several simplified procedures 
are being considered. At high rates of reactivity (large droplets at 
high pressures) , the combustion zone is located close to the droplet and 
the results of Ref. 3 can be used directly to discuss droplet response. 

At the other extreme, the analysis is considerably simplified when reaction 
rates are small and the droplet is simply evaporating. Finally, when the 
gas phase is quasisteady, consideration of response can be approached by 
using generalized transport relations which are also suitable for conducting 
droplet lifetime predictions. 

The activities in each of these phases of the investigation will be discussed 
in the following sections of the report. 

During this report period, several papers, theses and reports have been 
issued as a result of work done under this contract, as follows: 

C. B. Allison, "Burning Rate Response of Liquid Monopropellants to Imposed 
Pressure Oscillations," NASA CR— 134541, January 1974; also Ph.D. thesis. 

The Pennsylvania State University, University Park, PA, March 1974. 
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G. S. Canada, "High Pressure Combustion of Liquid Fuels," NASA CR-134540, 
January 1974, also Ph.D. Thesis, The Pennsylvania State University, 
University Park, PA, March 1974. 

G. S. Canada and G. M. Faeth, "Combustion of Liquid Fuels in a Flowing 
Combustion Gas Environment," Fifteenth Symposium (International) on 
Combustion, Tokyo, August 1974, also The Combustion Institute, Pittsburgh, 
PA, (in press). 

C. B. Allison and G. M. Faeth, "Open-Loop Response of a Burning Liquid 
Monopropellant," accepted for publication in AIAA J. 
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2* High Pressure Droplet Burning in Flowing Combustion Gases * 
2.1 Introduction 


The objectives of this phase of the investigation involve the determination 
of droplet vaporization rates and the characteristics (temperature, composition, 
etc.) of the liquid surface under forced convection conditions, at 
elevated pressures, in combustion chamber environments. 

Work during the present report period involved continued comparison 
of earlier theoretical results with measurements of droplet gasification 
in combustion gases at elevated pressures. Work was also begun on the 
development of an intermediate pressure model to allow more accurate prediction 
of the phase equilibrium at the surface of a burning droplet than is currently 
possible. 


2.2 Experimental Data Correlation 

Work reported in Ref. 1 involved comparison between theoretical 
and experimental droplet gasification rates and liquid surface temperatures 
over a range of conditions. The experimental portion of the investigation 
included measurements of surface temperatures and gasification rates for 
porous spheres supported in a flowing combustion gas stream leaving a 
burner. Data was obtained for the pressure range 1-40 atm. 

A number of theories were compared with these experimental results. 

The analysis involved the use of a spherically symmetic model for the gasi- 
fication rate in the absence of convection, in conjunction with a multipli- 
cative correction for the effect of convection, following a suggestion given 
in Ref. 4. The gas phase model of the combustion process involved a number 
of variable property treatments including; variable property— variable specific 
heat, variable property-constant specific heat and constant property versions* 
Phase equilibrium at the liquid surface was computed by convectional low 
pressure models (neglecting gas solubility and high pressure Poynting corrections) 
as well as a high pressure version for characteristics near the critical 
point, developed by Prausnitz and Chueh.^ 

The multiplicative correction has the form 

n/n o = 1 + 0.278 Re 1/2 Pr 1/3 (1 + 1.237 Re _1 Pr“ 4/3 ) " 1/2 

where 


n - convective burning rate 

n Q = burning rate with no convection 

Re = Reynolds number 
Pr = Prandtl number 


Following the suggestion of Ref. 4, the work in Ref. 1 employed approach 
conditions to evaluate properties in this expression. 
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Comparison between the various gas phase models indicated little 
difference in results as long as any assumed constant property was evaluated 
at an arithmetic average condition. The gasification rate predictions of 
the low and high pressure models were also essentially the same, with the 
exception that the high pressure model predicts a higher pressure for critical 
combustion. The overall results were not uniformly good, however, with errors 
ranging up to 50% between predicted and measured gasification rates. Agree- 
ment was generally poorest for high molecular weight hydrocarbons at low 
ambient oxidizer concentrations. 

During the present report period, this comparison was extended to 
include further measurements of droplet gasification rates in combustion gas 
environments (2) . Table 1 summarizes the range of test conditions for these 
two studies. 

In addition to the use of approach conditions in the convection 
correlation, a second method was considered based on average conditions in 
the flow field around a droplet. Eisenklam, et al. 6 have suggested a 
complex procedure for selecting average properties in a convection correlation. 
The method of Ref. 6, however, does not have correct asymptotic properties 
at low ambient oxidizer concentrations and it was not considered. For most 
fuels, the stoichiometry is such that the flame zone is located toward the 
outside of the flow field around the droplet. This is particularly the 
case at low ambient oxidizer concentrations. Therefore, based on this obser- 
vation, properties in the convection correction were selected at average 
conditions between the flame and the drop surfaces since this region comprises 
the bulk of the flow field. The average taken is based upon the average 
temperature and composition in the region. The average mixture properties 
are computed as described in Ref. 1. 

Some comparisons between predicted and measured gasification rates are 
shown in Figs. 1 and 2. Other examples may be found in Ref. 7. The rate of 
gasification has been normalized by the convection correction so that a range 
of test conditions can be shown on a single plot. Two sets of data are shown, 
one set normalized by the convection correction for approach conditions (dark- 
ened symbols) , the other set using average conditions (open symbols) . The 
theoretical curve is based upon a variable property-variable specific heat 
model of the combustion process in conjunction with both low and high pressure 
phase equilibrium models (1). Notably, the gasification rate prediction of 
these two models is essentially the same in the region in which they overlap 
(Fig. 2). 

In general, the use of average conditions in the convection correction 
results in a better prediction of the gasification rate. The improvement is 
greatest for high molecular weight hydrocarbons at low ambient oxygen con- 
centrations (lowest curve on Fig. 2). At higher ambient oxygen concentrations 
there is less difference between the two methods. This occurs since the 
approach temperature is a better approximation of the average of the droplet 
and flame temperature under these conditions. 

Measured and predicted burning rates are compared over the entire 
data range in Fig. 3. The data points shown involve the use of approach 
conditions in the convection correction. In order to avoid cluttering the 
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Table I 

Summary of Droplet Combustion Test Conditions 



Present 

- 2 

Faeth and Lazar 

Test Method 

porous sphere 

suspended droplet 

Ambient gas temperature (°K) 

600-1500 

1660-2530 

Ambient oxygen cone, (mole %) 

0-13.5 

0-41.8 

Pressure (atm) 

1-40 

1 

Diameter (cm) 

0.95 

0.11 

Approach Reynolds number 

40-680 

1.3-2. 2 



npAn/n 0 ) x iO^ (GM/SEC) 


6 
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600 


800 ' 1000 ' 1200 ' 1400 "“ISOO 
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Figure 1 Evaporation rate of methanol as a function of ambient temperature 
at atmospheric pressure. Approach Reynolds number range 46-174. 
Approach convection correction-darkened symbols; averaged convection 
correction-open symbols. 
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Figure 2 Gasification rates of n-heptane as a function of pressure for 
various ambient oxygen concentrations. Approach Reynolds 
number range 66-672. Approach convection correction-darkened 
symbols; averaged convection correction-open symbols. 
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Figure 3 Overall comparison of theoretical and experimental gasification 
rates in combustion gas environments • Lines A-A are data bounds 
when average properties are used in the convection correction. 
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figure, the data using average properties in the convection correction are 
not shown individually. The bounds using the data of Ref. l,with this 
latter procedure, are shown by the lines A-A. In the case of the measurements 
of Faeth and Lazar (2), the two procedures cannot be distinguished due to the 
low Reynolds numbers of these tests (Table 1) . 

Use of the averaged convection correction gave average differences between 
theory and experiment of ±30%, with maximum discrepancies ranging up to 70%. 

The use of approach conditions gives poorer results, particularly for the 
high Reynolds number conditions of the present tests. Finally, the present 
procedure for averaging is superior to that suggested in Ref. 6, since the 
method is continuous as the ambient oxidizer concentration approaches zero. 


2.3 Intermediate Pressure Model 


The low and high pressure phase equilibrium models have been found to 
represent surface conditions adequately when the ambient pressure is less 
than or greater than half the pressure required for critical combustion 
respectively. This is Illustrated in Figs. 4 and 5, taken from Refs. 8 and 9. 

A characteristic of the two models, illustrated in Figs. 4 and 5, is 
that they do not merge at low pressures. This occurs since the Prausnitz 
and Chuelr* model is adjusted to give a good representation of conditions 
near the critical mixing point and fails to present liquid properties 
adequately at lower pressures, where the conventional models are accurate. 

These results indicate the desirability of considering an intermediate 
pressure model that would provide a continuous transition between the two 
currently available versions. Since models in this region can be founded 
on a larger data base then the high pressure versions, they can provide more 
accurate predictions of high pressure phenomena. Since they formally merge 
with the low pressure versions, the Intermediate pressure models give 
information on the limit of usefulness for the simpler low pressure models. 
Under combustion chamber conditions, the loss of surface tension near the 
critical point results in shattering of the droplet, therefore, the inter- 
mediate pressure versions are likely to provide information over the bulk 
of the range where droplet combustion, via an envelope flame model, actually 
exists. Finally, the intermediate pressure models provide a good estimation 
of phase equilibria for complex hydrocarbon mixtures and thus have potentially 
greater usefulness for practical fuel types (as opposed to pure hydrocarbons) . 

The introduction of a new phase equilibrium model does not require a 
change in the theoretical treatment of the gas phase. In fact, the results 
to date indicate that the present gas phase model is adequate within the 
regimes pertinent to the low and high pressure phase equilibrium models 
(Fig. 2 for example). 

Of the various intermediate pressure models that could be used, the 
method outlined by Prausnitz, et al.^® has been selected for investigation. 
This procedure is well documented, and can be conveniently related to the 
current high pressure version. The detailed description of this model may 
be found in Ref. 10, in the following only a brief description will be given 
to indicate the main features of the method. 
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Figure 5 Liquid Surface Temperatures as a Function of Pressure for the Strand Combustion 
Experiments of Ref. 9. 
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For two phases to be in equilibrium, the total temperature and pressure 
must be equal in both phases, and the fugacity of each species in the system 
must be the same in both phases. 

v 

In the gas phase the fugacity, (J> p 

where p is the total pressure, Y* is the gas phase mole fraction of i and 
4>i is the fugacity coefficient of i. The model computes the fugacity co- 
efficient by the use of virial equation of state, taken to second order. The 
second virial coefficient is correlated in terms of individual and binary 
coefficients, which can be corrected for the presence of polar compounds in 
the system. 

In the liquid phase, the fugacity is given by f* = x^ 

where x-^ is the liquid phase mole fraction, y± is the activity coefficient and 
fj 1 is the fugacity in a reference state. The reference state for a condensible 
component (one that is ordinarily a liquid at the temperature of the mixture) 
is the fugacity of the pure liquid. If the mixture temperature is higher than 
the critical temperature of the component, Henry's constant is used as the 
reference fugacity, Wilson's equation is used to evaluate of activity coeffi- 
cient in both cases. This equation provides proper temperature dependence, and 
can be used for multicomponent mixtures while only requiring parameters obtainable 
from binary mixture data. The equation is also suitable for polar compounds. 

With correlations available for the gas and liquid phase fugacities, 
calorimetric quantities such as heats of vaporization can be determined as 
described in Ref. 1. With these results in hand, the phase equilibrium model 
and the gas phase model can be solved simultaneously to yield liquid surface 
temperature predictions for comparison with the experimental results. This 
requires physical constants for the pure species (which are available) , Wilson 
parameters, and Henry’s constants. The latter quantities are available in 
some cases, must be fitted from binary data (if this exists), or must be 
estimated. Reference 10 discusses various procedures for determining these 
quantities. 

At the present time, the computer program for the intermediate pressure 
model is being checked on the computer. The necessary property values for 
the correlation are being determined from the literature. When this work 
is completed, liquid surface temperatures will be predicted for the body of 
data available from earlier measurements (1, 2, 7, 8, 9, 11). The intermediate 
pressure model should overlap the range of validity of the high pressure model, 
allowing a check of the predictions of the two methods. 


2.4 Summary 

Work in the area of high pressure fuel droplet combustion has involved 
overall comparison of experimental and theoretical liquid surface temperatures 
and gasification rates in combustion gas environments. Predicted results are 
relatively independent of the gas phase model (when appropriate average properties 
are chosen). Using a convection correction with average properties between 
the droplet surface and the flame provides predictions within ±30% of measured 
values for typical combustion chamber conditions. 
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The low pressure phase equilibrium model is adequate for predicting liquid 
surface properties at total pressures less than half the pressure required for 
critical combustion* The high pressure model is more adequate at pressures 
approaching the critical combustion condition* 

Work has been initiated on the development of an intermediate pressure phase 
equilibrium model, to provide a continuous transition between the present low 
and high pressure versions. This work should be concluded in the next report 
period. 


3. Combustion Response of Monop gopellant Droplets 
3*1 Introduction 


The objective of this phase of the investigation involves the determination 
of the open-loop combustion response of monopropellant droplets to imposed 
pressure oscillations. Knowledge of the frequency ranges where the combustion 
response is maximum, allows the designer either to adjust combustion 
geometry to avoid having characteristic chamber frequencies in this region or 
provide damping tuned to this frequency range. 

The present work is a continuation of earlier efforts on obtaining a 
theoretical and experimental comparison for the response of a burning mono- 
propellant strand to imposed pressure oscillations* Using these findings as 
a basis, current activities involve examination of the more practical case of 
droplet combustion. By utilizing a model generally substantiated by strand 
combustion results, this step to the droplet combustion case can be approached 
with some confidence. 

Earlier work on the response of droplets has been confined to consideration 
of evaporation or bipropellant combustion. ^ Therefore, the present investigation 
should provide new information on the response characteristics of liquid fuels. 

In the case of liquid strand combustion, for a given fuel, the frequency 
response curve changes only as the mean pressure Is varied. ^ For monopropellant 
droplets, however, droplet size also influences the process introducing a new 
parameter to the problem. Going from large droplets to small ones, the 
reaction zone becomes more influenced by curvature, and lies at proportionately 
greater distances from the surface for the small droplets. This implies that 
large droplets (where the reaction zone is close to the surface) behave somewhat 
like strands , while small droplets behave somewhat like evaporating droplets 
without reaction. Therefore, the present work provides information in the 
branch region limited on both sides by earlier results. 

Initial work on this portion of the investigation involved, first of all, 
direct consideration of the large droplet limit in the context of the earlier 
strand combustion model. Then, using parameters determined in Ref. 3, the 
transient analysis is recaste into a form appropriate for droplet combustion. 

The above approach is limited to a frequency regime where the characteristic 
time of oscillation Is small in comparison to the droplet lifetime. At low 
frequencies, a different procedure must be used, more analagous to the 
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vaporization response calculations based upon droplet lifetime (12, 14, 15). 

At this limit, the gas phase is essentially quasisteady providing a substantial 
simplification of the analysis. This approach can also be linked to the 
more complete perturbation analysis at the limit of a quasisteady gas phase. 
Therefore, the analysis treats these cases in a unified manner. 


3.2 Large Droplet Limit 

A burning liquid strand, itself, is analagous to a segment of an infinitely 
large liquid droplet. This analogy persists for a range of droplet sizes, allowing 
direct use of the results of Ref. 3 in this "large droplet" regime. 


The actual size range of the "large droplet" regime will be considered 
later. In general terms, this is a condition where the reaction zone is close 
to the surface. For droplets of this type, ignition occurs before the bulk 
liquid has heated to any degree, and there are temperature gradients in the liquid 
phase. gi nce the reaction zone is close to the surface, the gas phase 
combustion process is largely uninfluenced by forced convection (unless very 
massive convection is present) .13,16 p n this situation, the combustion process 
is completely analagous to liquid strand combustion. 


As an example of typical response characteristics. Fig. 6 shows the response 
for a hydrazine strand. 3 The response is given as the real part of the acoustic 
admittance as a function of the dimensionless frequency of the imposed pressure 
oscillations. This parameter is equivalent to the response parameters defined 
in Ref. 12 for liquid fuel combustion, at the large droplet limit. 

With increasing frequency, two positive peaks are observed in the acoustic 
admittance plot. The first peak occurs near the characteristic frequency of the 
thermal wave in the liquid phase, and is a result of the interaction between , 
the combustion process and liquid phase transient effects. The second peak 
occurs at frequencies near the characteristic frequency of the gas phase thermal 
wave, and is a result of interaction between the combustion process and gas 
phase transient effects. 


Combustion instability is usually associated with values of the real part 
of the acoustic admittance greater than a value on the order of unity . ^ Table 2 
summarizes the fequency ranges where this criteria is satisfied for the results 
of Fig. 6. The results have been converted to dimensional form with the two 
frequency ranges corresponding to the liquid and gas phase transient regimes. 

At atmospheric pressure, only the gas phase transient range has sufficient 
amplification to drive instability, while at pressures greater than 10 atm., both 
regimes are of significance. 


Typical frequency ranges for combustion instability and typical sizes for 
large droplets must be considered in order to relate the findings of Table 2 to 
practical combustion chambers. Combustion Instability is usually associated with 
the frequency range 500-30,000 Hz.^-^ At low pressures, Table 2, indicates 
that this range would involve gas phase transient effects; while for pressures 
greater than 10 atm. interaction of the combustion process with the liquid phase 
thermal wave is the more significant phenomena. These findings suggest that 
while the use of a quasisteady gas phase assumption at low pressures might be 
adequate for mean burning characteristics of monopropellants, potentially important 
contributions to the combustion response can be overlooked with this approach. 



REAL PART OF 
THE ACOUSTIC ADMITTANCE 



Figure 6 Computed Real Part of the Acoustic Admittance as a Function of Frequency and 
Pressure for Hydrazine • 
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Table 2 

Critical Frequency Ranges for Hydrazine 
Strand Combustion 


Mean Pressure Liquid Transient Gas Transient 

(atm) Range (HZ) Range (HZ) 


1 - 300-1500 

10 approx. 15 4500-21500 


100 


450-1650 


>225,000 
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The droplet size appropriate to the "large droplet" limit can be examined 
by estimating of the position of the reaction zone. Spalding* s theory of mono- 
propellant droplet combustion (as presented in Ref. 19) provides a convenient 
approach to this problem when strand combustion rates are available. From this 
theory, it can be shown that the dimensionless flame position is given by 
(see Ref. 19 for notation) 

a 2 = (r /r ) 2 = [ 2 + r + (r 2 + 4 r ) 1/2 ]/ 2 r 

c c e L J 


where 

p v d c /r = (1 + 0 ) X In [ 1 + (1- 0 )/6 ] 

eeep e c e 

e = t /t 

e e c 

6 = L/c T 

P c 

For hydrazine, at pressures greater than atmospheric pressure, the liquid 
regression rate is given by^ 

v = 0.01254 p 
e 

where v g is in cm/sec and p is in atmospheres. 

Using property values for hydrazine given in Ref, 3, the droplet size can 
be computed for various flame positions and pressures. These results are given 
in Table 3 for dimensionless flame positions of 1.1 and 2.0. 

In Table 3, for a =1.1 the decomposition flame is only 10% of the droplet 
radius away from the liquid surface. Droplets larger than the values listed in 
Table 3 for this a should behave as "large droplets" (this estimation is 
preliminary, however, pending the results of the next phase of the investigation). 
Significantly, at the higher pressures, the "large droplet" regime covers the 
range of droplet sizes usually encountered in rocket engine applications. This 
indicates that the magnitude of the response (Fig. 6), the frequency range 
(Table 2) and the droplet size (Table 3) are all commensurate with the 
requirements for combustion instability under rocket engine conditions. 

In Table 3, for cr c = 2, the flame is located relatively far from the 
droplet and there is only weak interaction between the gasification of the liquid 
and the decomposition flame. Judging from earlier analyses of droplet evaporation 
and combustion response under these conditions, *- 2 the liquid transient response 
is likely to be much reduced at these conditions. The transient gas phase 
response, however, could still contribute at these conditions. Therefore, at low 
pressures, the unsteady gas phase may still be an important factor. 

The drop size range between the two limits shown in Table 3 is a branch 
region where the liquid phase response should be increasing toward the large 
droplet limit. This clearly comprises a significant technological range of 
droplet sizes, indicating the importance of considering actual droplet response 
as opposed to just the large droplet limit. This, coupled with the possible 
importance of gas phase transient effects for small droplets at low pressures. 
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Table 3 

Hydrazine Droplet Sizes as a Function of Flame 
Position and Pressure 



1 5250 290 

10 525 29.0 


100 


52.5 


2.9 
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provides motivation for the more detailed analysis of droplet response 
considered in the next section. 


3.3 Droplet Response Analysis 

3.3.1 Description and Assumptions of the Model 

The previous work has suggested the need for specific examination of the 
response of burning droplets. Since response is sought, as opposed to droplet 
lifetime calculations, the possible importance of gas phase transient effects 
must be considered. The previous work on strand combustion response has been 
substantiated by experiment, therefore, the present theoretical effort to 
extend the results to droplets can be approached with some confidence. 

The configuration under consideration involves a droplet that is burning 
in the absence of convection, yielding a spherically symmetric flow field. 

While convection is important under actual combustion chamber conditions, 
there is evidence to indicate that the region where the droplet has a low (or 
zero) relative velocity with respect to the ambient gas is a very critical zone 
in determining combustion instability characteristics.^ For simple evaporation 
without decomposition, the response due to velocity effects has been found to 
be quite low. 20 in addition, the decomposition process near the droplet sur- 
face reduces the influence of convection over a fairly wide range of conditions. 1 51 
This occurs since convection only influences the process when the outer edge 
of the flow field interacts with the reaction zone, a situation that is limited 
to very weak reactivity or very high Reynolds numbers. Therefore, neglecting 
convection for monopropellant droplets puts fewer limitations on the practi- 
cality of the calculations, than would be the case for bipropellants. 

3 20 

Similar to earlier response studies, * the ambient pressure is assumed 
to be oscillating with a wavelength that is long in comparison to the 
dimensions of the combustion field of the droplet. The period of oscillation, 
however, is assumed to be short in comparison to the total droplet lifetime so 
that consideration does not have to be given to large changes in the position 
of the surface during a pressure oscillation. This assumption allows the 
analysis to proceed while only examining oscillatory solutions. The low 
frequency regime where this breaks down will be discussed later (this low 
frequency regime is largely associated with the range of frequencies where both 
liquid and gas phases are quasisteady) . 

Since the combustion rate of monopropellant droplets is much higher than 
bipropellants, the usual assumption of a constant liquid phase temperature 
(in the zero order) is less valid than for bipropellants. Therefore, the 
presence of mean liquid phase temperature gradients must be considered in the 
analysis. Examination of constant mean temperatures can be made, in order 
to include the conventional steady combustion model, by equating the bulk 
liquid temperature to the wet bulb temperature at a given pressure. 

A sketch of the theoretical model is shown in Fig. 7. Notation may be 
found in Table 4. The process is examined at an instant of time when the 
droplet radius is r g . Formally, this radius is taken to be fixed so that 
the mass flux of fuel is time varying. This actually corresponds to porous 
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Symbol 



n 

P 

q 

qs 

r 

R 

t 

T 


w 

Y, 


a 

3 

Y 

r 

6 

5 i 

e 

6 


TABLE 4 

Notation for Open-Loop Response Analysis 
Description 

vapor pressure parameter, Eq. (21) 
parameter, Eq. (11) 
preexponential factor 
specific heat 

binary diffusivity 
activation energy 
heat of formation 

n 

heat of vaporization, Eq. (21) 
vapor pressure parameter, Eq. (21) 

mass flow per unit solid angle 
molecular weight 
reaction order, Eq. (10) 
pressure 

energy of reaction 

liquid phase heat flux, Eq. (58) 

radial distance 

gas constant 

time 

temperature 
radial velocity 
reaction rate 

mass fraction of species i 

thermal diffusivity 
parameter, Eq. (21) 
specific heat ratio 
parameter, Eq. (55) 
reaction parameter, Eq. (10) 
parameter, Eq. (4) 

amplitude of oscillatory pressure, Eq. (29) 
combined variable, Eq. (23) 
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X 


t 


Vi 

P 


9 


tl 


a 


* 


0) 

U) 


s 


thermal conductivity 
stoichiometric parameters, Eq. (5) 

parameter, Eq. (9) 

density 

parameter, Eq. (69) 
parameter, Eq. (66) 
parameter, Eq. (52) 
frequency 

frequency for strand case, Eq. (53) 


Subscripts 


f 

F 

i 

s 

0 

1 


liquid 

fuel 

center of droplet 
liquid surface 
steady state quantities 
first order quantities 
ambient conditions 


Superscripts 

* dimensional quantities 
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sphere combustion, however, the two cases are equivalent as long as the 
density in the liquid phase is large compared to that of the gas phase and 
the period of oscillation is small in comparison to the lifetime of the droplet. 
When this is true, the motion of the surface with respect to the mean surface 
position is negligible and can be ignored. Exceptions to this assumption 
arise near the critical point and the present analysis is not valid in this 
regime. (Several other aspects of the model would also have to be modified 
in order to consider near critical conditions) . The response portion of the 
analysis is invalid at frequencies having an oscillation period comparable to 
the lifetime of the droplet. 

For generality, gas phase transient effects are included in the model, 
since results discussed earlier indicate that could be important at low pressures. 
The effect of variable properties are also included, so that the model is 
equivalent to the earlier strand combustion analysis (3). 

The remaining assumptions of the analysis are similar to those of Ref. 3. 

They are as follows: 

1. The process is spherically symmetric with a Mach number much less 
than unity and negligible body forces. Inertial and viscous terms in the 
momentum equation are neglected. 

2. The reaction process is premixed and laminar. A one-step, 
irreversible chemical reaction takes place in the gas phase and any time 
lags associated with the chemical reaction itself are negligible, i.e. the 
chemical reaction is locally quasisteady and obeys an Arrhenius equation valid 
under steady state conditions. Chemical reaction is neglected in the liquid phase. 

3. Radiation heat transfer is neglected. 

4. The gas phase is taken to be an ideal gas and the Lewis number is assumed 
to be unity. 

5. All gas diffusion coefficients are equal, all molecular weights are 
equal, all gas phase specific heats are equal and constant, the gas phase thermal 
conductivity is independent of composition and varies linearly with 
temperature, and the liquid is composed of a single chemical species having 
constant properties. 

6. The combustion products are assumed to be insoluble in the liquid phase 
and the gas phase fuel mass fraction at the liquid surface is given by the 
Clausius - Clapeyron equation. As in the case of the gas phase reaction, the 
equilibrium at the surface is assumed to occur rapidly in comparison to other 
transient effects in the system. 

7. The wavelength of any periodic pressure disturbance is assumed to be 
long compared with the diameter of the zone involving active combustion. Con- 
sideration of the momentum equation, along with Assumption 1, then implies 
that pressure is only a function of time. 

A discussion of the applicability of these assumptions is provided in 
Ref . 3 . 
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3.3.2 Governing Equations 


The notation for the following analysis is provided in Table 4. 
to simplify notation, define the following dimensionless variables: 


In order 


r = r*/r* 
s 


t - t* X* /P *C *r * 2 
°°o o°o p s 


v *0 *c *r */X * 

r «o p s °°o 


© = m*C */r *X * 
p s °°o 


( 1 ) 


T - T*/T w o * p = p*/p Q * p 

X « X*/X * D = D*/D * 

CDQ OOQ 




* 


where the subscript °° 0 refers to a zero order (steady state) quantity, evaluated 
at infinity. The conservation equations are now considered in the liquid and gas 
phase, followed by the boundary conditions. 


Liquid Phase r < 1 

The equation of conservation of mass is 


p f r V fr = m f 


( 2 ) 


where m f is the mass flow rate per unit solid angle and is only a function of 
time. Conservation of energy is 



3T 
fr 3r 



3r 


, 2 3T. 

(r 3F> 


where <5^ is the ratio of thermal diffusivities 


(3) 


5 


f 




( 4 ) 


Gas Phase r > 1 


The stoichiometry of the one-step gas phase reaction is assumed to have the 
following form: 


N N 

yui- ^EZ V [1] (5) 

i=i 1 i=i 1 
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Overall conservation of mass in the gas phase is 


9£ 

3t 


+ 


2 

r 


9r 


(pr 2 r r ) 


0 


( 6 ) 


It is also convenient to define a gas phase mass flow rate per unit solid 
angle, analagous to Eq. 2, as follows: 


m 


2 

pr v 


In this case, m is a function of both radial distance and time. 
Conservation of species becomes 
3Y. 


P( 


+ v 

9 t r 


3 V 1 9.2 4 , 

~ 1 • ~2 9 ? (r PC - 5 F> + V i w - i - 1 . N 


m 


( 8 ) 


where 


N ± = <V - V }/ V 

w = A exp (-E/T) 


, C *B*r * 2 T 
A = _g s °°o 

X * 

°°o 


r p 



M*R*T * 
L- ooq _ 


n 

3 


and 


(9) 

(10) 

(ID 


E = E*/R*T * 

OOQ 


Finally, conservation of energy is 




9T. 1 9 

3r ; r 2 9r 


, 2 .3T. 
(r ^ 


+ 0 


Y - 1 


■> # + 


dt 


q w 


where 


q ■ 



V.'J/CV ' C *T *) 
i F p °°o 


( 12 ) 


(13) 


(14) 


is the dimensionless heat of reaction, and p is only a function of time. 
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The assumption of an ideal gas having a unity Lewis number and with thermal 
conductivity propertional to temperature; provides the following relationship 
between properties: 

T = X = pD = p/p (15) 


Boundary Conditions 

Since the liquid is incompressible, it has a constant temperature at the 
origin 


r - 0 T - T (16) 

In the case of steady combustion at the wet bulb temperature of the droplet, 

Ti would be equal to the wet bulb temperature. However, since monopropellants 
have higher gasification rates than most bipropellants, cases exist where 
temperature gradients are present in the liquid phase throughout the lifetime 
of the droplet. Therefore, T i is taken as a parameter for the analysis in order 
to investigate the effect of these gradients on the response. 

Aside from continuity of temperature and pressure at the surface of the 
droplet, four other conditions must be satisfied. These may be summarized as follows: 


r - 1 


P £ v rf 


Pv 


r 


X 


f 


3T 
3r 

s- 



s+ 


pD 




p v r [(1-3) T + L] 
P \ (Y p - 1 ) 


Y f = ± exp (-L f /T) 

where 


(17) 

(18) 

(19) 

( 20 ) 


3 = c */c * 

f p 

L = L*/R*T * 

a = a*/p * 

*o 

L t - L f*/ RT »o* 


( 21 ) 
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Equation (17) represents conservation of mass at the surface. Equation 
(18) is the conservation of energy relationship, where the heat flux reaching 
the surface goes into the heat flux in the liquid phase and the heat of 
vaporization of the gasifying fuel. Equation (IS) follows from the 
insolubility of the product species in the liquid phase which implies that 
the total mass flux of the non— fuel species must be zero at the surface of 
the droplet. Equation (20) is the Clausius- Clapeyron equation relating the 
fuel mass fraction to temperature at the surface. Formally, Lf and L should 
be equal, however, the vapor pressure curve of the fuel is often better 
represented by a value of Lf slightly different than that of the heat of 
vaporization so that provision has been made for this possibility. 

Far f ran the droplet, all the fuel is consumed and the temperature is independent 
of distance 

' r = co Y = 0 T = T (t) ' (22) 

F 

where T(t) is a known function depending upon the specified pressure variation. 

3.3.3 Schvab-Zeldovich Transformation 

The solution can be simplified by introducing a Schvab-Zeldovich variable 
based on the fuel mass fraction 


0= q Y_ + T 


(23) 


In terms of this variable, Eq. (8) is replaced by 


, 30 . 30 v 1 3 n 30 % + ( Y “ h *E 

p< 9t + v r 37 > ^ 3? (X 37> + 


dt 


(24) 


Gas phase conservation of energy is unchanged, however, the reaction rate 
becomes 


w = A q~ n p n T S_n ( 0 - T) n exp (-E/T) (25) 

The boundary conditions at the liquid surface, Eqs. (19) - (20) are as follows: 

r = 1, 

* 5 < e ~ T - ) -) = pv (6 - T - q) (26) 

3r s+ r 

0 _ t = 51 (-L /T) (27) 

P t 

The boundary conditions far from the droplet become 

r = CO 

0 = T » T (t) (28) 


The remaining equations are unchanged. 
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3.3.4 Perturbation Analysis 

Following the approach used in Ref. 3, a linear perturbation analysis 
is conducted. It is assumed that the pressure oscillates with small amplitude, 
e, and the dependent variables are represented as follows: 

p (t) ® 1 + £ e iwt 

T(r,t) * T Q (r) +eT 1 (r)e il,),: 

0(r,t) = 0 o (r) + £ 0 1 (r)e iWt (29) 

m(r,t) = m Q (r) +ei^{r)e^ 10t 

. , . . . . ioit 

m f (t) = ra fo + e m fl e 


Substituting Eqs. (29) into Eqs. (1-28), the resulting equations are separated 
into like powers of ee^. Equation (15) is also employed in order to reduce 
the equations to forms involving only 0 and T. 

The zero order, steady state, equations are as follows: 

m, = m = const (30) 

to o 

In the liquid phase, r ^ 1 


9 dT “fo 

(r 2 -£> - P f S f 
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The first order equations are as follows. In the liquid phase. 


d , 2 Ui l. 

3F (r cGT> ", 


= const 

dT - *£0 

p £ « £ dr 


. 2 A dT 

iar T _ fl o 

6 f X 1 " p f 5 f dr 


In the gas phase, r > 1 


2 d 1 2 2 

T ", iiirr T. = - iwr T 

o dr 1 o 


H 5 dG l d0 T ? d0 

I f(tl "~r“) - m T - r lufl, = (m T- + m.T ) ~ 

o dr o dr o o dr 1 o 1 1 o dr 

, 0 d0 , „ d0 0 

- T 4— (r 2 T -^) - T n ~ (r 2 T -p- 2 -) - r 2 ( Y - — -) 1 w T 

o dr v 1 dr 1 dr v o dr Y o 


dT dT dT 

T %- (r T -p^) - m T -r^ - r dioT, = (m T. + nLT ) 
oar o dr o o dr 1 o 1 1 o dr 

, , dT , . dT „ ,, , 

-T f (r 2 T -J—2-) - T 4— (r 2 T —•) - r 2 -) i to T 

o dr 1 dr 1 dr o dr Y o 

T 

- r 2 q 1_n A T o 1+6 ‘ n (0 q - T Q ) n exp (-E/T q ) [n+(l+6-n)^ 

U J 0 

0, - T ET- 

+ n (— “t + ~2 ] 

o o T 

o 


The boundary conditions on these equations are 


r = 0 


T 1 = 0 


r - i 

dT dT dT 

x £ *r>_ ■ T „ &rK . + T 1 sr* .. - 4 o (1 - B) T 1 - “i 


s- 


s+ 


s+ 


[(1-3) T o + L] 

d(9 - T ) d(6 - T ) 

T o dr >s+ + T 1 V" >s+ - “o (9 1 " V + A 1 

L T 

0 1 - T x = a q exp {-L f /T Q } [ - 1 ] 


r ^ 1 

(39) 

(40) 

(41) 

(42) 


(43) 

(44) 

(45) 

(46) 

(47) 


2 


(48) 
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From Eq. (28), 9, and T. must approach a constant value from the surface of the 
droplet. Solving Eq. (53) for large r it can be shown that 


r = co 



(49) 


which is the form for isentropic flow. 

The solution of Eq. (31) can be obtained immediately for the boundary conditions 

T = T. at r = 0 and T = T at r ■ 1, as follows: 
o i o os 


T - T. 
o i 


T - T 
os i 


exp 


m 


fo 


p f 5 f 



(50) 


Equation (50) provides an expression for the steady state temperature distribution 
in the liquid phase. 

In the gas phase, Eq. (32) can also be integrated once without difficulty. 
Applying the boundary conditions at the droplet surface, Eq. (35) and (36), and 
determining the liquid phase temperature gradient from Eq. (50), yields 


2 

r T 


d^ 

dr" 


o 


oo 


(51) 


where 


\p =0 + (A/p<5 ) (T - T) + L - B T - q (52) 

r o o K j os i os 

In general, the solution of the remaining equations must be obtained 
numerically* For the zero order problem, Eqs. (33) and (51) must be integrated. 

There remain four boundary conditions that have not yet been used, among Eqs. (35)- 
(38). Three of these provide boundary conditions on the differential equations, the 
fourth determines the eigenvalue m^. 

The first order equations, Eqs. (39) — (43) » are linear, requiring at total of seven 
boundary conditions. Equations (44) -(49), plus the comp at ability conditions 
at the liquid surface, T^ - T lg , provides eight boundary conditions. The extra 
boundary condition determines the eigenvalue mf^, which can be related to the 
response. 

The solution of the entire system of equations, Eqs. (29-(52), represents 
a formic&ble task, even on the computer. Numerical calculations along this line 
are currently in progress and the results will be reported later. Under certain 
circumstances, however, the computations can be simplified considerably. The 
formulation under these conditions will be undertaken in the following. 

3.4 Quasisteady Gas Phase Analysis 

3.4.1 General Analysis 

Consideration of the strand combustion case in Ref. 3, indicated that transient 
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gas phase phenomena contributed to the combustion response when the strand 
dimensionless frequency, w was of order unity. The present dimensionless 
frequency is related to the earlier value as follows. 


co = 


m *C *r 


<»o 


* 2 

-) 


U) 

s 


(53) 


At high rates of reaction, the factor in parenthesis in Eq. (53) is a 
constant. This quantity can be evaluated from the hydrazine droplet combustion 
data of Ref. 21. The results indicate that at the limit of high reactivities, 
important gas phase transient effects should not be encountered for dimensionless 
frequencies less than 1000. In this region, the gas phase may be assumed to be 
quasisteady with all important transient effects confined to the liquid phase. 

The quasisteady regime also encompasses the low frequency region where the 
time of oscillation is an appreciable fraction of the lifetime of the droplet 
and the previous perturbation analysis is invalid. Therefore, In eve op ng 
relations for the quasisteady limit, an objective is to obtain a uni re 
treatment both by perturbation analysis and eventually by droplet lifet e 
calculations. 

In order to achieve these objectives, the gas phase is analyzed, neglecting 
time derivatives. The boundary conditions of the problem are the droplet 
surface temperature and the ambient conditions. With this information, the 
mass fraction of fuel at the surface is also known through the Clausius- 
Clapeyron equation. The desired results of the analysis are the heat and mass 
fluxes at the droplet surface. This information can then be used to determine 
both response and lifetime characteristics of the droplet. 

The assumptions of the analysis are the same as in the earlier treatment, 
except for additionally neglecting gas phase transient effects. Under these 
assumptions, the conservation equations are as follows. 


Af-ir -o 


(r 2 X — - m T) + r 2 q w 
dr dr 


= 0 


where 


r - 6 + q o + 1 -B T - q 

a o 


^ dT. 
m dr 


s- 


(55) 

(56) 

(57) 

(58) 


The boundary conditions on these equations are 


T = T 

s 


r = 1 


(59) 
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0 = 0 - 3S, exp (-L,/T ) + T 

s p f s s 


( 60 ) 


X 


-) 

dr'' 


s+ 


= m 


[q s + (l-e)T s + L) 


r = so 


(61) 


T = 6 - 1 (62) 

The definitions of 0, w, etc. are the same as before and Eq. (15) provides 
a representation of the physical properties. The mass flux per unit solid 
angle, m, is a constant throughout the gas phase. The quantity q , defined 
in Eq. (58) is the heat flux to the liquid phase after the heat of vaporization 
requirements at the surface have been met. 


Equations (55) and (56) must be integrated numerically in general. Given 
T g and p, Eqs. (60) yields 0 S . The equations only require one additional boundary 
condition, Eqs. (61) and (62) provide three. Therefore, the additional boundary 
conditions determine the eigenvalues A and q s . In this manner, the solution 
of the equations determines the functions 

m = m (p, T g , A) (63) 

q s = q s (P, T s , A) (64) 


where A carries the Influence of droplet size. 


3.4.2 Gas Phase Analysis for Small A 

For very small droplets, the parameter A is small, and the effect of 
reaction can be neglected. At this limit, the droplet simply evaporates. 

The solution of the gas phase equations can be completed analytically under 
this conditions, as described in the following. 

With A equal to zero, Eq. (56) may be integrated once, yielding the 
following after applying the boundary condition of Eq. (61). 

(65) 

where 

* = T + q + L - B 1 (66) 

s s 

Eliminating spatial derivatives between Eqs. (55) and (65) and then integrating 
yields an expression for ® , as follows: 

6 + q +L-BT - q T + q +L-BT 
S S _ 2$ s 

0 e + q„ + L-8 T - q T+q+L-gT 
s’s s n s n s s 

after applying the boundary condition at the droplet surface. 


( 67 ) 
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Equation (65) is also readily integrated. After applying the surface 
boundary conditions the following is obtained: 


m (1 - — ) = T - T - a In [(T+a)/(T + a)] 

r s s 


( 68 ) 


where 


a = q + L - B T 
s s 


(69) 


The eigenvalues, m and q are determined by applying the outer boundary conditions, 
Eqs. (62), to Eqs. (67) and (68). This yields 


m = (1-T ) - a In [ (1+cr)/ (T + a)] 

S s 

i s - 1 ' (1 ' e T s + l) 

s s 

where it is recalled that 

8 S - f ‘‘W + T s 


(70) 

(71) 


(60) 


Equations (70), (71) and (60) provide the necessary transport parameters at the 
droplet surface. They represent analytical forms of Eqs. (63) and (64) at the 
limit of A equal to zero . 


3.4.3 Perturbation Analysis 

If the liquid phase equations are perturbed, Eq. (50) still represents 
the zero order solution. The steady state surface termperature is then 
related to the bulk liquid temperature, as follows: 


X r dT , 

«s. - r sr > - W ( (T os - T i> 

o s- t 


(72) 


Therefore, given T ± and p, Eqs. (64) and (72) can be solved to yield T og . The 
zero order gasification rate, m Q » can then be obtained from Eq. (63). 

The first order liquid phase equations are also unchanged from Eqs. (39) and 
(40). The boundary conditions, for a given value of A, are 


r = 0 T 1 = 0 


(73) 
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fl 
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X dT l 

f j — ) = m q , + m. q 

<*r o si 1 so 
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where 


♦ 9m \ _ . 9m. 

a i = *r> „ T is + 

S p T s 


(75) 
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T i + 
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*. 

9p 


(77) 


The derivatives in Eqs. (76) and (77) are known framthe general gas phase solution, 
evaluated at T , p . Solution of Eq. (40), with the boundary conditions of Eqs. 
(73)-(77), allows tRe eigenvalue m^,to be determined. 


When analytical forms are available, as is the case when A * 0 (Eqs. 70 
and 71), the derivatives in Eqs. (76) and (77) can be determined exactly. For 
finite values of A, only numerical results are available for the functions 
represented by Eqs. (63) and (64). In this case, the derivatives can be 
determined from plotted numerical results using p and T g alternately as parameters. 
This procedure essentially corresponds to the Zeldovich method (22), which has 
recently received some attention for the calculation of solid propellant 
combustion response (23-25). Unlike the solid propellant case, however, the 
present calculations must include droplet size as a parameter, through A. 

At the present time, the quasisteady gas phase analysis has been programmed 
on the computer for the general case of any value of A. Calculations are also 
proceeding at the limit of small A. These numerical results will be completed 
in the next report period. 


3 . 5 Summary 

During this report period, results for burning liquid monopropellant strand 
have been applied directly to determine preliminary droplet response characteristics 
for large droplets. The generalized analysis has also been completed allowing 
for transient phenomena in both gas and liquid phases. 

The system has also been analyzed at the limit of a quasisteady gas phase. 

In this case, simple analytical forms have been obtained, as well, at the limit 
of an evaporating droplet. 

At the present time, calculations are in progress to determine the open loop 
response function for monopropellant droplets (12) . Initially work is considering 
response under the limitations of a quasisteady gas phase. The transient gas phase 
response will follow after the simpler case has been considered. The results of 
these calculations will be presented In the next report period. 
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